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ABSTRACT 

A formalism for jet noise prediction is derived that includes the refractive 'cone of 
silence' and other effects; outside the cone it approximates the simple Lighthill format. 

A key step is deferral of the simplifying assumption of uniform density in the dominant 
'source' term. The result is conversion to a convectedwave equation retaining the basic 
Lighthill source term. The main effect is to amend the Lighthill solution to allow for 
refraction by mean flow gradients, achieved via a frequency-dependent directional factor. 
A general formula for power spectral density emitted from unit volume is developed as 
the Lighthill-based value multiplied by a squared normalized' Green's function (the 
directional factor), referred to a stationary point source. The convective motion of the 
sources, with its powerful amplifying effect, also directional, is already accounted for in 
the Lighthill format: wave convection and source convection are decoupled . The 
normalized Green's function appears to be near unity outside the refraction dominated 
'cone of silence', this validates our long term practice of using Lighthill-based approaches 
outside the cone, with extension inside via the Green's function. The function is obtained 
either experimentally (injected 'point' source) or numerically (computational 
aeroacoustics). Approximation by unity seems adequate except near the cone and except 
when there are shrouding jets: in that case the difference from unity quantifies the 
shielding effect. Further extension yields dipole and monopole source terms (cf. Morfey, 
Mani, and others)) when the mean flow possesses density gradients (e.g., hot jets). 

INTRODUCTION 

Lighthill, in his seminal papers 1 - 2 posed the problem of flow noise in terms of a wave 
equation for a 'virtual medium at rest': the actual flow was effectively incorporated into 
right hand side terms; these were interpreted as sources of sound. Although the equation 
is exact, approximations to these source terms have the effect of suppressing sound 
convection (hence refraction and shielding) by the mean jet flow 3 * 4 . (Some effects of 
refraction were pointed out by Powell^ even before this connection was made.) 
Equivalent equations for a moving medium have been put forward 6 ' 10 ; they allow for 
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the sound convection. Of these, Lilley's equation 1 ** has received much attention: it has 
been developed by Mani 11 ' 13 , Balsa 14 - 1 * 1 , and others 16 ' 22 into a quantitative predictive 
formalism for properties of jet noise. It entails a formidable derivation or calculation of a 
Green’s function for highly idealized models of a jet flow. By contrast, the Lighthill 
procedure, as developed by Ribner 23 , Pao and Lowson 24 , and others based on 
Csanady 7 - 25 " 27 , is relatively simple. These various formalisms appear to yield 
comparable predictive accuracy outside the 'cone of silence' opening downstream of the 
j et 28,29 Within this cone the Lighthill-based theory fails completely— it predicts no 
attenuation— whereas the Lilley-based theory exhibits good to poor accuracy, depending 
on frequency. 

It is now well known 3 - 4 - 7 that expansion of the basic Lighthill source term leads to extra 
terms that may be shifted to yield a convected wave equation: it was the implicit 
discarding of these extra terms that was cited above. In 1977 the expansion was 
exploited 28 - 21 * to demonstrate a considerable equivalence between the Lighthill- and 
Lilley-based approaches outside the 'cone of silence'. In the course of the present study it 
was realized that the. dominant residual source terms, as finally modelled, were fully 
equivalent to the Lighthill term. Taken together, these findings led to the notion that the 
Lighthill solution could he amended to allow for convected wave effects: extended in 
simple fashion into the 'cone of silence'. The basis of the extension lies in replacing the 
ordinary oscillatory Green's function, e 11<r /4jtr, by that for the convected wave equation. 
These are both for a stationary point source, in contradistinction to the Lilley procedures; 
this decouples convection of the sources and the sound waves, permitting the cited 
simplifications in the theory. The following analysis develops the mathematical 
implementation. Since the Green's function is frequency-dependent, the formalism is 
directed toward the power spectral density, in particular, from unit volume. Then 
quadratures can provide the full jet spectrum and the broadband noise intensity. All of 
these are direction-dependent. 

The survey paper, reference 29, from which the present study evolved, illuminates certain 
facets only briefly touched on herein; moreover, it displays graphically several 
comparisons of theory and experiment merely cited here. For a fuller understanding and 
perspective as to how the present notions relate to other theories of jet noise, that paper 
should be consulted as well. 
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PARTIAL OVERVIEW 


Lighthill manipulated the conservation equations of fluid dynamics into the form of a 
wave equation forced by nonlinear terms in the unsteady velocity v; on the right hand 
side; these were interpreted as sources of sound. For practical prediction, only the 
dominant source term was retained, to yield 


1 d 2 p V 2 p _ YL PJirJ 


(^pvjvj 

CoW " " 3xi9xj (la> 

As a further approximation, the fluid density was approximated as uniform (p = p 0 ) in 
the source term, resulting in 

JL <^2 y2 n _ n 3 2 VjVi 
2 0 t 2 -Vp - Po 3 0 
c 0 J 


(1) 


The early assumption of uniform density is premature, however: it completely suppresses 
refraction of the sound by the flow 3 ’ 4 . Herein, gradients in the density figure in the 
expansion of the original source term of (la), after which the density is replaced by its 
mean. This restores the refractive capability in terms of a convective wave equation: in 
the case of a transversely sheared mean flow it takes the simple form 

■&XI (2a) 


^ • * 2 P - PoS* 


(where U2ac is defined by 


Pu2ac _ dp 
P° Dt 3x2 


(2b)) 


as a refinement of the unconvected wave equation (1): the degree of approximation in the 
step p=p 0 here is . less. This bears a close relationship to Lilley's well-known 
equation 10 ’ 12 * 13 . It is, however, vastly simpler, with the special feature of being forced 
by the same source term as in Lighthill's uniform density format, if, compatibly with p = 


Po. 3vj/3xj is taken to be zero. The convection implicit in D/Dt is by the jet velocity, 
approximated locally as U(x 2 ); this is the time average of the instantaneous velocity vj, 
the fluctuations being designated Uj. 


The solution of Eq.(2) for a single frequency component of the pressure p is obtained via 
a Green's function, the right hand side being replaced by a stationary oscillatory point 
'source'. Guided by an elegant relation due to Balsa 15 (but like Lighthill, he used a 
moving point source), we develop a general expression for the power spectral density of 
the sound pressure emitted from unit volume of the flow as 
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dO(xlco) 
d Vol 


«= IGI 2 * 4-dimensional Fourier transform of space-time correlation of 


(p 0 d 2 ViVj/dxjdxj) (3a) 

where IGI is the absolute value of the oscillatory Green's function, and it applies to a large 
distance x (far field). (An earlier attempt by Schubert 8 * 9 , in retrospect oversimplified, 
may be mentioned.) Eq.(3a) can be reduced to the much neater form 

^ Vot^ °e k 4 * IGI 2 * 4-dimensional Fourier transform of space-time correlation of 


(poV K - 2 ) (3b) 

where K is a wave vector (with absolute value k) governing the phase of G, and v K is the 
instantaneous velocity component in the direction of k. This reduction parallels a similar 


one in a space-time domain formalism (for p 2 ) due to Proudman 45 . 


The four dimensional transform may he reduced to a one dimensional transform (time => 
frequency); this involves the widely used approximation of the turbulent ’eddies' as 
'acoustically compact'. This is a great simplification, but at a significant price in 
accuracy: the amplifying effect 1 * 2 of eddy convection is progressively overestimated with 
increasing jet velocity3» 4 »30,38 


An alternative version of Eq.(3) is obtained by transformation of the correlation to a 
convected frame of reference, following Chu^U 1 . This resembles equations of Lilley 111 
and Balsa 1 3 based on a moving point Green's function in connection with a moving 
frame. (Cf. also references 24, 38, 46-48 for developments in terms of a four dimensional 
Fourier transform.) In either case, use of the moving frame simplifies extraction of the 
amplifying effect of convection of the source pattern. 


For the simple wave equation (1 ), the magnitude of the Green's function at a large 
distance x is simply (1/4tix). For the convected wave equation (2), on the other hand, IGI 
will have a form approximating ( 1/4tcx) times a directional factor 8 * 9 . In terms of these 
results, a solution ofEq.(I) corresponding to (3b) is obtained by replacing tc^lGI 2 in 
Eq.(3b) by k 4 (4rcx)' 2 , where k = co/c. Then, on taking the ratio of the two equations, 

d<t>(x|co)/dVol = [(K 2 /k 2 ) (4;tx)IGI] 2 dd>(x|o))/dVol (4) 

convected wave eq. Lighthill wave eq 
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where the factor in brackets is defined as IG^yl and for the purposes herein is designated 
the normalized Green's function. 

Eq.(4) generalizes the power spectral density from unit volume based on the simple 
Lighthill wave equation to account, via a more realistic wave equation, for the influence 
of flow. It exhibits the result as the Lighthill-based value multiplied by a squared 
'normalized' Green's function. The Green's function multiplier incorporates the flow- 
acoustic interaction: it yields a frequency-dependent alteration of the directional pattern. 
The interaction has been described as a refractive effect 32 ’ 33 due to mean flow or sound 
speed gradients and also as a shielding effect 1 L 12 * 15 due to the mean flow itself. The 
term 'Lighthill-based value' in the above relation is not restricted to a formulation in 
terms of 4-dimensional transforms. All valid estimation schemes used in the past are 
encompassed; these are usually based on an approach in the space-time domain, rather 
than the wave number-frequency domain. 

Equations (2), (3), and especially their corollary Eq.(4 ) are the key results of this paper. 
They provide the basis for reinterpretation of early procedures of the author’s group — 
experimental 32 ’ 33 and numerical 8 - 9 ’ 34 — that effectively evaluated an approximation to 
the Green's function. The normalized magnitude, IGa'I, appears near unity outside the 
refraction dominated 'cone of silence': this validates our practice of using Lighthill-based 
approaches outside the cone, with extension inside via the Green's function. 

A further extension in an Appendix yields dipole and monopole source terms (cf. 
Morfey 35 , Mani 13 , Michalke and Michel 49 - 5() ) when the mean flow possesses density 
gradients (e.g., hot' jets). 


MODIFIED LIGHTHILL EQUATION 
Unconvected Wave Equation: Virtual Medium at Rest 

Lighthill's wave equation 1 - 2 is equivalent to 

1 8 2 p y2 _dWj 1 02 p d 2 p 
c 0 2 dt 2 P “ 3xjdx ; + c 0 2 dt 2 ’ 3 t 2 


(5) 


wherein the pressure has replaced the density that he used as the dependent variable; 
additionally, the fluid has been approximated as inviscid and non-heatconducting. The 


5 



right hand side is interpreted as a spatial distribution of sources of sound. In the usual 
approximation, the last two terms are taken to cancel, and the fluid density is taken to be 
a constant (p = p G ) in the first. This leads to 


1 3 2 p 

C 0 2 3t 2 


^ - Poli* 


( 6 ) 


A rich literature 4 - 29 has dealt with applications of this equation for the prediction of 
properties of jet noise. 


The early replacement of p by p 0 to yield (6) is premature, however; it has the effect of 
suppressing wave convection (and refraction) by the flow 3 ' 4 . This is shown most simply 
by expansion of the original first source term under the specification of a unidirectional, 
transversely sheared, mean flow U(X 2 ). (This is known as the 'locally parallel’ 
approximation when applied to a real spreading jet. That more complicated case is dealt 
with in Appendix A, along with the derivation of additional source terms 13 ' 35 that arise 
in, e.g., heated jets, from density gradients. ) The instantaneous local velocity is written 
as the mean plus a perturbation uj, 

vi = Ui + Ui ; Uj = (U(x 2 ),0,0), (7) 


and the expansion changes (5) to 

1 3 2 p v2 _ ^P 11 ^ ? 3U3pU2 1 3 2 p P 2 p 

c () 2 3t 2 ’ P ~~ dxjdxj + “3 x 2 0X1 + Co 2 3t 2 ’ Dt 2 

using the definition D/Dt = 3/3t + U3/3x] 

as a convective derivative following the mean flow. Both (5) and (8) are exact. 


( 8 ) 

(9) 


Convected Wave Equation: Actual Medium with Flow 

At this point we approximate D 2 p/Dt 2 as c' 2 D 2 p/Dt 2 \ where c is a local time-average 
sound speed (Appendix B). On shifting the term to the left hand side, equation (8) goes 
over to 


1 D 2 p 

— -)Dt 2 
c l 


- V 2 p 


3 2 pujUj 3u3pu2 
3xj3xj + ^3x 2 3xi 


( 10 ) 
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Thus in place of the unconvected wave equation (6), we now have the more accurate 
convected wave equation ( 10). The difference is traceable to the deferral of the 
incompressible flow assumption in the first source term of Eq.(5): this had suppressed the 
expansion of 9 2 p/9t 2 into the convective form D 2 p/Dt 2 . 


This is not yet in final form. The second source term, involving mean flow shear 9U/9x2, 
is linear in U 2 . For this reason, it has been argued 8 ' 10 ’ 36 that the term participates in wave 
propagation and so should be on the left hand side. (See Goldstein 37 , pp. 389-391, for a 
further discussion.) This applies, however, only to a small acoustic (or compressible) 
component associated with wave propagation. Within a subsonic flow the 
overwhelming part of uj is induced by the turbulence vorticity; being small compared 
with the soundspeed, it may be approximated as incompressible. Thus we split off the 
acoustic component and place it on the left hand side. (The acoustic component of the 
first source term, on the other hand, is of higher order and may be neglected). Thus it is 
at this point that we may justifiably apply Lighthill's approximation p= p Q to the 
remaining right hand side source terms. Consistently, the uj are taken to have zero 
divergence. Equation (10) then goes over to 

^ • V’P - PO^ * 2PoIS 

c z J 

where U2ac is defined via the momentum equation 


CD 


Pu2ac _ 3j>_ 

Po Dt “ ‘5x2 

The corresponding equation for the Green's function, G(x,ylo>), is 


icot 


(11a) 


(12) 


and the solution of the pair of equations (11a) and (12) for p is Ge' icot . 


Recall now that Eqs. (8) to (12) relate to turbulence uj superposed on a transversely 
sheared mean flow Uj = (U(x 2 ), 0, 0), Eq. (7). For this scenario we examine the Lighthill 
source term p 0 9viVj/9xjXj of Eq.(6), where vj = Uj + uj. On carrying out the 
differentiation, the term expands exactly into the two terms on the right hand side of 
Eq.(l 1). Thus that equation is exactly equivalent to 



auau^ 

ru dx2 


r-fj 92viVj 

p - Po 9xj9xj 


(13) 
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In the expansion we have used the incompressibility relation dvj/Oxj = 0 implied by taking 
p = p 0 ; as noted, the approximation of incompressibility is made only in the source terms. 

This is a major new result, approximated as Eq (2) in the OVERVIEW. It is slightly 
generalized here, with c 2 (x) replacing c 0 2 . (A further generalization in Appendix A, 

(i) replaces the restrictive Eq. (7) by a more realistic mean flow U(x) on the left hand 
side, and (ii) allows for mean flow density gradients; these give rise to important 
additional source terms in the cases of hot jets and jets of nonambient gas 13 - 35 * 49 - 50 .) 

Now according to Lighthill’s arguments 1 ’ 2 , p 0 dvjvj/5xjxj is a valid source term for flow 
noise. From the above it is clear that both the right hand source terms of Eq.(l 1), being 
equivalent for the sheared mean flow, are likewise valid source terms for that scenario: 
the ’shear noise’ term 2p o 0U/dx2)0U2ac/dxi) as well as the 'self noise’ term p o 0ujuj/dxixj. 
Furthermore, the velocities uj therein are incompressible to the same approximation as in 
Lighthill’s term. All of this supports the earlier argument concerning the split of the shear 
term in Eq.(l 1): a compressible acoustic component has been excised (before the step 
p=>p G ) and placed on the left hand side, leaving an incompressible component on the 
right hand side (presumably much larger). The shear term with U2ac on the left hand side 
serves in a propagation role, whereas the shear term with U2 on the right hand side serves 
in a generation role. 

Equation (11) with deletion of the last term, together with Eq. (11a), are roughly 
equivalent to Lilley's equation 10 - 12 ’ 13 (he combined them into a single third order wave 
equation 37 ). However, Lilley, in effect, moved the entire shear term to the left hand side 
wave operator. The 'shear noise' source term so lost from the right hand side is of major 
importance. This was shown indirectly by Ribner 4 - 23 , and later more directly by Pao and 
Lowson 24 in terms of Lighthill's equation. Aside from the arguments of the last 
paragraph, credibility is afforded by comparisons with experiment: the 'self noise' 
spectrum, downshifted in frequency, shows a predicted match to the 'shear noise' 
spectrum (Fig.l). We will return to this point later. 

The successive modifications of the wave equation for flow noise, proceeding from 
Eq.(5) (virtual medium at rest) to Eq.(13) (actual medium with flow) may be summarized 
at this point. Deferring the approximation p=Po, as in the derivation of Eq.(13), leads to 
propagative terms that may be moved to the left hand side. This is a conversion into a 
convective wave operator, making explicit the role of the mean flow in convecting the 
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sound waves. On the other hand, the convective role of the mean flow is only implicit in 
Eq.(5): it resides in the gradients of p in the source term. These gradients are suppressed 
in the approximation p=p 0 prematurely applied in Eq. (6). At very high frequencies wave 
convection should reduce 7 ’ 1 1.54,55 the we u known amplification predicted!- 2 , in the 
absence of flow, for moving sources; the reduction is referred to as 'fluid shielding'* *• 14 . 
But at jet noise frequencies the most dramatic effect has to do with the flow gradients: the 
radial decrease in mean-flow velocity in the jet; this turns or refracts the sound rays away 
from the axis to yield the well known 'cone of [relative] silence’. 

The allowance for refraction in the generalized Lighthill equation (13) comes at a price: 
the point source solution (Green's function, utilizing eq. (1 1) as well) is enormously more 
complex than that for the simple wave operator of Lighthill's equation (1). It has been 
found, however 8 - 9 , that there is near agreement at the larger angles from the jet axis. This 
implies that the simple Lighthill solution for rms sound pressure will approximate that 
predicted via eq. (13) outside the 'cone of silence'. The basis for this is elaborated 
below. 


FORMULAS FOR POWER SPECTRAL DENSITY 
General Relations 

The formulas that follow are based on the Green's function for a stationary, oscillatory 
point source. The approach parallels that of Balsa! 5 based on a moving, oscillatory point 
source, with missing steps being inferred. We seek the power spectral density O(xlto) of 
the radiated sound pressure dictated by equations (1) or (13); they may be written 
symbolically as 

L[d/d t, d/dx; a(x)| p(x, t) = Q(x, t) (14) 

where L may be either the unconvected wave operator of Eq (1) or the convected wave 
operator of Eq (13). The a(x) are the coefficients; for the convected wave operator they 
allow for the local mean flow, taken as U(x 2 ), and a space-variable sound speed. (In a 
generalized version of Eq. (13) in Appendix A, U(x 2 ) goes over to U(x),and p 0 to p(x), 

where p is a local time average. Equation (14) applies to this version as well.) 
Correspondingly, 

L[-io), d/Bx; a] p(xlto) = £>(xlco) (15) 

where p and Q are defined in Fourier transform pairs: 
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/>(xlco) = 2^jr Jp(x,t) e ia)l dt; p(x,t) = Jp(xlco) e' ,(0t da) 


-oo 
oo 


-oo 

oo 


Q(x|co) = J Q(x,t) e ia)t dt ; Q(x,t) = Jg(xlG)) e' ia)t dt 

-oo -oo 

Equation (15) has a solution 

OO 

p(xko) = J G(x,yl(o) 0(ylw) d 3 y 


(16) 

(17) 


(18) 


in terms of a Green's function G(x,ylco) that is the solution of 

L[-ico, d/3x; a] G(x,yloo) = 8(x - y) (19) 


The (two-sided) power spectral density of the sound pressure is evaluated as 


<J>(xlco) = (p(xlco) p*(xlG))> 


( 20 ) 


where ( ) signifies an ensemble average. Inserting (18), with y replaced by y' and y", 

OO 

OO f 

<D(xlco) = I J G(x,y'ko)G*(x,y’’l(o)<0(y'i(o) 0*(y"la))>d 3 y’ d 3 y" (21) 

-OO - 00 

oo 
OO - 

= j jG(x,y+^/2|a))G*(x,y-^/2|(o)(0(y+^/2|a))0*(y-^/2|a)))d 3 ^d 3 i; (22) 


where £ = y' - y" and y = (y' + y")/2. 


Reduction for Far Field 

The Green's function of Eq. (22) may be written in the form 

G(x,y 1(0) = |G(x,y ko)| e'M^y to) (23) 


We now restrict x to the far field defined by Ixl being very much greater than both lyl and 
the largest wavelengths of concern; y is limited to the region of nonzero source strength 
Q. A sufficient approximation for the phase, which seems to be implied in Balsa's 
Appendix 15 , is then 
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\|/(x,y'lco) = \j/(x,y"l(0) - K»(y'-y") (24) 

where the wave vector 

K = — (^y Y)far field (25) 

is proportional to CD. There seems to be the further reasonable assumption, which we 

make also, that the variation in amplitude of G is negligible compared with that of the 
phase as £ of eq. (22) ranges within the source region Q. Then ^ may be dropped in 
comparison with y in the amplitude so that 

G(x,y+^2|co) G*(x,y-^2|co) - |G(x,yl<B)| 2 e‘ i,c ^ (26) 

The other factor in (22) is the frequency-domain correlation 

K(y, Slco) = (G(y+^/2|co)0*(y-£/2|co)> (27) 

This is the Fourier transform of the time-domain correlation 

R(y* S.t) = (Q(y+^/2, t+x) Q(y-^/2, t)) (28) 

(which is independent of t); specifically, 

OO 

R(y , m = ^ J e ia)X R(y, ^,x)dx (29) 


Insertion of equations (26) and (29) into (22) yields 

OO OO 

OO 

<D(xlco) = i J IG(x,y|(D)l 2 d 3 y J J R(y£,x) e - iK ^ +i(OT d 3 £dx (30) 

-OO -00-00 

Equation (30) is the desired general result for the power spectral density &(x\ co) (cf. 
Eq.(3) of the Overview). The inner integral can be recognized as a four-dimensional 
Fourier transform of the two-point space-time correlation R(y, £,x). Alternatively, it is a 
three-dimensional transform of the cross-spectral density R(y, £|co). This transform 
multiplied by the square of the amplitude of the Green’s function (frequency domain) has 
a simple interpretation: it is the contribution of unit volume of the sources Q to the power 
spectral density of the sound pressure radiated to the Field point x. 
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Moving Reference Frame 

Experimentally the correlation R(y,^t) in a jet flow has a form describing a moving, 
fluctuating pattern. This is dealt with most neatly by transforming to a reference frame 
moving with the pattern convection velocity, taken as U c . (But the Green's function, 
unlike that in Balsa's relation 15 , still refers to a point source at rest.) Following Chu 3l > 
we take 

= ^ - UfT ; U c = (Uc,0,0) (31) 

so that - k*U c t (32) 


and reexpress R in terms of using (3 1 ), as 


Rm(jS^mT) - R(y^.t) 

Then in Eq. (30) 

oo oo oo oo 

J jR(y4,T)e-i^+i««d34dT = J J R,„(y£m.*) 

- 00-00 - 00-00 


(33) 

(34) 


since the Jacobian of the transformation is unity. We may further define 


G3 = 0) - k*U c 


(35) 


as the effective source frequency in the moving frame to yield an observer frequency (0 
at x in the stationary frame (far field). Inserting these last two equations converts Eq. 
(30) into 

OO oo oo 

4>(xIcd) = 2 n JlG(x,yl 0 ))l 2 d 3 y J J R m (y£ m /0 e iK *S in +iroT d 3 £ m dx (36) 


This is an alternative form for the power spectral density; it is more useful in that the 
effects of source convection are more easily brought out. In Eq. (30) the space-time 
source field correlation is referred to a stationary coordinate frame and is designated R. 
In Eq. (36 ) this same correlation is referred by transformation to a coordinate frame 
moving with velocity U Cx _ and is designated R m . In both cases the two points being 
correlated are stationary. 
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Source Term of Form a(y)8 2 bjj/8yjdyj 

Suppose the (monopole) source strength distribution has the form of the right hand side 
of Eq. (13), slightly generalized: 

Q(y, t) = a(y)3 2 bjj(y,t)/5yidyj (summed over i, j = 1 to 3) (37) 

It can be shown that the Fourier transform of this is 

fi(ylco) = a(y)a 2 ^(yl(0)/ay i 3y J The 

The cross spectral density of the source field between points y' and y” is then 

dyjdyj dy koy | 

which is equivalent to 

<fi(y'lo»j2*(y"la))) = a 2 (y) (38) 

where % = £i, ^3 = y' - y", if a(y')a(y") is approximated a 2 ((y'+y")/2) = a 2 (y). 

Insertion of equations (26) and (38) into (22) yields the power spectral density of the 
radiated sound pressure as 


<t>(xlco) = 


IG(x,y|co)l 2 a 2 (y) - ~~~ — e-‘ K ^ d 3 £d 3 y (39) 


which is replaceable by 

oo 

O(xlco) = J IG(x,y|co)l 2 a 2 (y) KiKjK K K|(/?y(y , la))6* / t/(y"la))> e' iK *^ d 3 ^d 3 y (40) 

-oo 

But the source cross-spectral density designated by ( ) is the Fourier transform of a two- 
point correlation: 

OO 

(bij(y'\(i))b*kl(y”\(o)) =^~ J<bij(y’, t+x)b kl (y", t)) e^ dx 

271-00 

so that the power spectral density of the sound generated by sources of type (37) may be 
written 

OO 

O(xl( 0 ) = | IG(x,y|co)l 2 a 2 (y) KjKjKkK|(bjj(y\ t+x)bki(y", t)> e-' K *^ +ia}t dxd 3 2|d 3 y (41) 


Suppose we identify a(y)bij(y,t ) with the Lighthill source term pfyjv/vy (wherein the 

uniform p G is generalized to p(y) to allow for a space-variable time average density: cf. 
Appendix A). Then, on taking KjKjKkKi inside the ( ) of Eq.(41), we have terms like 
KiKjbij = KiKjV|Vj . But the summation KjVj is k times the component of v along k, 
which we designate v K . Writing k as the magnitude of k, this is summarized as 
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( 42 ) 


KjKjbjj = KjKjVjVj = K 2 V K 2 

It follows that, in Eq.(41) the summation 
KiKjK k Ki(bij(y', t+T)b k i(y”, t)> = k 4 v K - 2 ((y\ t+T)v K 2 (y", t)> = k 4 R K (y£,T) (43) 
where R K is a two-point correlation of v K 2 . With this replacement Eq.(41) simplifies to 

oo 

r - 00 

d>(xlco) = ^ J »c 4 p 2 (y)IC(x,ylo))l 2 d^y J | R*(yfrz) e -'^+'«w d^dx (44) 

-OO -00-00 

where, it is noted, the two points being correlated in R K are referred to a stationary 
reference frame. Correspondingly, Eq. (36) becomes 

00 

00 

1 f °° 

<J>(xlto) = 2 ^ J K 4 p 2 (y)IG(x,y|co)|2 d 3 y j J RKm(y4m^)e- iK ^m +iG5T d 3 ^ m dx (45) 

-OO -00-00 

where the correlation R K - m designates the same correlation with respect to the moving 
frame (velocity U c ); it is obtained from R k - via the transformation (31). Eq.(45 ) is a key 

result stated in words as Eq.(3b) of the Overview (hut generalized with p(y) replacing 
Pc <)■ 

These latest results for the power spectral density O(xlco) may be put in perspective: they 
are all expressed in terms of the Green's function for a stationary, oscillatory point source 
in an arbitrary flow U c (x). Equations (30) and (36) refer to a general source strength 
function Q; Eqs.(43) and (45), on the other hand, refer to a source strength of the 

generalized Lighthill form Q = p(y)<D 2 vjVj/0yjdyj (summed), to which p(y)K 2 v k 2 is 
equivalent. 

It is noteworthy, as he pointed out, that this double divergence form of Q implies that the 
sources are of quadrupole nature. In the time domain this was associated with an 
operator d 2 /dt 2 in the far field format; in the present wave-number domain the factor k 2 
plays an equivalent role. 


14 



Virtual Medium at Rest 

This is Lighthill's scenario 1 - 2 : the fluid flow is incorporated into the source term. The 
original region of flow is now treated as a 'virtual medium at rest'. For this case the 
oscillatory Green’s function is simply 

G(x,ylto) = e*'*-yl ; k = 0)/c o (46) 

When the observer point x is in the far field (cf. after Eq. (23)), a sufficient 
approximation is, with x = Ixl, 

G(x,ylco) = 4 ^ eikx ‘ ikx#y/x; IG(x,ylto)|= ^ (47) 

Thus, in this case, the vector K = -(V y \j/)r i , r field of Eq. (25) may be identified with the 
wave vector k given by 


k = k = kx/x = o)x/c 0 x ; k = lk| (48) 

Also, with the sources convected parallel to the xi-axis, U c = (U c ,0,0), Eq. (35) yields 

OJ = 0)|l -(U c /c 0 )(x/x)l = co(l-M c cosG) s to© (49) 

This is just the Doppler shifted source frequency that will yield an observed frequency 0) 
at x. 

Reduction to Chu's Equation 

Equation (45)— correlation referred to moving reference frame, Q of form 
Po3 2 viVj/3xj9xj-may be applied to Chu's scenario 30 - 31 by specializing to a fluid at rest. 
This implies invoking Eqs. (47) to (49). His power spectral density d>i(x|co) is one-sided 
(limited to positive values of (0), so that it is twice our d>(xko). He notes further that, 
since Rkm(y^m.'t) is even in ^ and X, the Fourier transforms may be replaced by cosine 
transforms. Equation (45) then takes the form, after doubling, 

oo oo 

d>l(xl(o) = P - 9 - — JlG(x,ylo))l 2 d 3 y Jcos k»^ m d 3 £ m jRkmfy^m.'Ocos G5T dx (50) 

-oo -oo -oo 

But we note that, from the equation for the time delay t m * above Eq. (8) of ref.31, we 
have the equivalence 
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O)0Xin* - k*^ m (51) 

employing our Eq. (49). Inserting this and the appropriate values of IGI, and G3 from 
Eqs. (47) and (49) 

co oo oo 

^(x!®) = 7 7 - 3 - 4 ~ 2 j d3 - v Jcos (O0T m *d 3 ^m jR km (y,^ m .t)cos co0tdx (52) 

1 OTC Cq X -00 -00 -00 

This recovers Chu's result for the one-sided power spectral density (Eq.(2.2.7) in ref.(30) 
and Eq.(8) in ref. 31). 

The simplifying assumption that the turbulent 'eddies' are 'acoustically compact' permits 
a marked simplification: cos (O0T m * in Eq. (52) (and similarly e _,K *^ m in Eq. (45)) can be 
replaced by unity. The effect is to reduce four dimensional Fourier transforms to one 
dimensional transforms: the space-wavenumber transform becomes merely a volume 
integral. The replacement is explained along with other observations concerning Eq. (52) 
in the following quotation from Chu3 () (annotations shown in brackets; our notation and 
reference numbering): 

"Firstly, we can identify 03 as the frequency in the turbulence; the corresponding radiated 
frequency co is then the Doppler-shifted frequency (i.e., (0 = 03 /0). This is a logical 
result which one could have obtained on physical grounds.. Secondly, Lighthill's 
criterion for neglecting retarded time shows up automatically in the cos co0T m * term. 
According to Lighthill, retarded time can be neglected if (0L/c o is small so that the eddy 
size L is small compared with the wave-length of the sound that it generates 
l 'acoustically compact'! . If this condition is met. ..then for ^ m < L the term cos G)0T m * 
can be approximated as unity. Thirdly, if retarded time [difference across an 'eddy'l is 
neglected, Lighthill's convection factor ( 1-M C cos0)~ 5 , which accounts for the main 
effect of convection at limited speeds is exhibited as a vertical shift (1-M C cos0)~ 4 plus a 

Doppler shift (1-M C cos0)~' in the power spectrum This concept is a low-speed version 

of a similar idea presented in ref. 4. Fourthly, although the Lighthill's convection factor 
is not applicable for high-speed convection because of its singularity where 1-M C cos0 = 
0, |an example given by Chu, discussed above Eq. (60)] will show that Eq. (52) is still 
valid for high-speed convection if retarded time [difference] is not neglected. In fact, this 
moving frame integral with proper account of retarded time posses a zereo that exactly 
cancels the (1-M C cos0)~ 5 singularity and replaces it by a nonsingular convection factor 
(cf. also refs. 3 and 38)." 
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Actual Medium vs Virtual Medium 

Here we compare the noise emission from the actual medium, allowing for the effect of 
the fluid flow on propagation, to that predicted for the virtual medium at rest. We will 
show how the former differs from the latter (the Lighthill scenario) in being an extension 
to allow for flow-acoustic interaction effects: e.g., refraction that bends sound rays away 
from the jet axis to create a 'cone of silence' opening downstream. And we will show a 
close to asymptotic approach to the Lighthill case outside the cone of silence. 

It will be convenient to restrict attention to the power spectral density emitted from unit 
volume at y. For the scenario of Eq. (45) ('actual medium'), with the appropriate Green's 
function and anticipating that k may be approximated as k at jet noise frequencies 
(Appendix C), 

a3< ^ ICd) = Po 2 k 4 IG(x,y|w)n2{^- J } Rkm(y,UA)e- ik ^m +iC3T d3Udx} (53) 

dy J 2 h-oo_oo 

The corresponding result, using the respective Green's function and wave vector k = k for 
a virtual medium at rest (Lighthill format), is 

OO °° 

[ 830(xltd) ]vM = p Q 2 k 4( 1/47tx) 2(J- j j R km (y,UA)e- ik ^ m +i05x d 3 Udt} (54) 
ay 5 2 n -oo.oo 

The ratio of these two can be put in the form 

a3< ^ (X ' C ° ' = [(4nx) IG(x,y|co)l| 2 [^^]vm (55) 

d y-’ dy-’ 

It will simplify discussion if we refer to the factor in brackets as a 'normalized Green's 
function’ 


IG/v(x/x,ylco)l = (4nx)IG(x,y|co)l (56) 

so that (54) may be written 

- = IG/v(x/x,yl(0)l 2 F (57) 

dy- 5 dy J 

with IGwl dependent solely on the direction of x by virtue of the 4ttx normalization and 
the 1/x decay of IG I in the far field. 
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Normalized Green's Function, |G,v|: Single Choice for Entire Jet 

The 'normalized Green's function' designated in Eq. (56) is nominally a function of source 
location. In this section we develop evidence that the value for a single choice of source 
position y on the axis, designated y 0 , may serve as an effective average. 

This was concluded in the context of a series of experiments 32 ’ 33 on the far-field 
directivity pattern of a 'point source' of sound immersed in a subsonic jet. Except for 
some uncertainty as to the accuracy of simulation of a point source, the measurements 
effectively yielded values of/Gp normalized by the value at 90°. This result was deemed 

equivalent to IGyvl 2 , on the ground that the effect of wave convection was expected to be 
nil at 90°; that is, IGI should reduce to (l/4jtx) there. It was found that the geometric 
average of the t)— dependence for symmetric off-axis positions +t> or -0 differed little 
from that of the on-axis source position for a given yi. Further, the variation with yi was 
small. This justifies referring IG/vl 2 to a single location, which will greatly simplify both 
utility and interpretation. 

Noting that the direction x/x may be designated t),<p in polar coordinates for a general, 
non-round, jet, the effective average of the squared 'normalized Green’s function' may be 
defined as IGyvl 2 ; thus 

IG jV y rc |iO))l 2 = ((4rcx) 2 IG(x,y|w)l 2 > ave (58) 

where y rc j is a representative value of y i along the jet axis. (The dependence on cp, of 
course, disappears for a round jet.) According to the arguments above, this refers to an 
average over y. In practice, however, it would be used as a surrogate for a weighted 
average; that is, equivalent to a single average value of IG(x,y|co)l 2 taken outside the 
integral of Eq. (45) to replace the value inside that varies with y. The use of a surrogate 
single Green’s function has also been the practice in solutions of the Lilley equation 12 ' 13 ; 
Mani 13 referred to the same experiments cited above as justification. Replacement of 
the y-dependent Green's function of Eq. (57) by the single y-independent Green's function 
of Eq. (58) allows immediate integration of (56) into 

cb(xlco) = IG/v (d,(p> y<)lw)l 2 !<I>(xIco)]vm (59) 

This states that the same frequency-dependent /G,v/ 2 that applies to unit volume is, to a 
sufficient approximation, applicable to the jet as a whole. 

DIRECTIVITY OF JET NOISE 
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|Gjvp Defines Flow-Acoustic Interaction. 

The form ofEq. (59) tells us that the flow-acoustic interaction effects, refraction and 
shielding, are embodied in the normalized Green's function /Gn/ 2 . This factor modifies 
the smooth directional pattern of intensity otherwise predicted (IGyyl 2 taken as unity). 

Both experimental (injected 'point' source) 32 ’ 33 and numerical (computational 
aeroacoustics) 8 ’ 9 ’ 34 attempts to evaluatelGwl 2 display this. The most striking effect is a 
progressive reduction of intensity within a 'cone of silence’ or 'refraction valley’ opening 
downstream along the jet axis: this is owing to the sound having been refracted outward 
by the jet velocity gradients. For filtered jet noise there is a similar 'cone of silence' that 
is matched by the experimental \G^\ 2 impressively well 28,29, at each frequency (Fig. 2). 

For the computed IG^I 2 there is a qualitative match 8 ’ 9 , with fairly good agreement at low 
frequencies and Mach numbers. However, the computed depth in decibels of the 
'refraction valley' is much exaggerated at high frequencies and high subsonic Mach 
numbers. We would attribute this overprediction to the use of the mean flow only in the 
convective terms of the wave operator (the turbulence being averaged out): the neglected 
large scale, low frequency, velocity field distortions would spoil the assumed 
axisymmetry during the transit time of a wave. It is intuitively evident that perfect 
axisymmetry is required to yield a very high refractive attenuation along the axis (e.g., 

54 decibels in Schubert's 8 - 9 most extreme case.) 

The frequencies for which geometric acoustics is applicable are many-fold higher 8 ’ 9 than 
those of jet noise. It was found 8 * 34 that at these high frequencies the computed valley 
depth is grossly exaggerated (e.g., 90 dB prediction at M=0.3). Despite this, many 
studies (e.g., references 7 - 17 ) attempt to quantify the 'cone of silence' via geometric 
acoustics. 

For heated or cooled jets, or jets of foreign gases, sound-speed gradients come into 
play 8 - 9 - 32 " 34 . Heating enhances the outward refraction, hence increases the depth of the 
refraction valley. Cooling has an opposite effect. If cold enough, the temperature 
gradients could dictate inward refraction strong enough to overpower the outward 
refraction imposed by the velocity gradients. This would give rise to some ’focusing' 
enhancement of noise intensity along the jet axis: IGnI 2 should exhibit an axial lobe in 
place of a 'cone of silence'. This expectation was dramatically confirmed in the 
experiments of Grande 29 - 33 : the enhancement lobe was 9 dB at 3000 Hz for an M=.l 12 
jet of nitrogen at - 1 80 C. An almost identical lobe was found in his measurements of the 
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jet noise in a narrow filter band tit the same frequency (Fig.2). Schubert’s approximate 
numerical calculations of IGa'I- showed a similar, albeit exaggerated, lobe 8 * 9 . 

This similarity of behavior between jet noise and experimental and computational 
approximations to IGa - I 2 , for both ambient and cold jets, shows that IG^I 2 dominates the 
intensity pattern in a conical downstream zone about the jet axis; elsewhere it is near 
unity. These facts, together with Eqs. (57) and (59), permit the following interpretation: 
|Gn|, the normalized Green's function, serves to extend the Lighthill-based calculations 
into the refractive zone near the jet axis: the 'cone of silence', or focused lobe, as the 
case may be. Approximation of IG^I 2 by unity outside this zone yields the Lighthill- 
based pattern: this seems an adequate approximation except when there are shrouding 
jets. 

The calculations need not be formulated in terms of the four-dimensional Fourier 
transform of Eq. (54), despite its figuring in the derivation: other formulations deemed 
to be equivalent may be used. In practice, relatively crude approximations— implicitly 
for the correlation R K - m — have been used. Some approximate developments of this kind, 
normally resulting from an analysis in the space-time domain, are discussed in a later 
section and in Appendix D. 

Convective Amplification and Basic Directivity 

Lighthill 1 showed that convection of the sources (but not the sound waves) gives rise to 
amplification of the noise in the downstream direction. To quantify this, Chu 30 modeled 
Rkm as it Gaussian function in the convected frame and carried out the four dimensional 
Fourier transform to evaluate Eq. (52) or the equivalent Eq. (54). The resulting power 
spectral density per unit volume was then integrated over frequency to yield the broad 
band radiated intensity per unit volume. This displayed the M c — > 0 intensity as 

multiplied by a convective amplification factor, C 5 , given by 

C* 5 s [(1- M c cosr» 2 + a 2 M c 2 ]' 5 / 2 (60) 

in terms of the source pattern convection Mach number M c . This factor was originally 
obtained by Ffowcs Williams 3 ** and by Ribner 3 via different formalisms; it is a 
refinement of the Lighthill-based relation 39 (1- M t cost^)* 5 that is nonsingular when 
M c cos f) = 1. (The basis for this was elaborated under Eq. (52).) The directional pattern 
of intensity radiated by the sources (the Basic Directivity), before multiplication by this 
factor, was tacitly taken to be spherical. 
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Ribner 23 made allowance for the two source terms of Eq.(l 1), leading to a nonspherical 
basic directional pattern; he chose a more defensible (but still somewhat unrealistic) form 
for Rkm- it was modelled as compatible with isotropic turbulence superposed on a 
transversely sheared mean flow. This could be expressed as a 'narrow band' directional 
spectrum 

C-5 H C[a(Cto)+ b(Cco)(cos 4 fl + cos 2 fl)/21 } 

= C' s * Basic Directivity (quasi-elliptic) (61) 

emitted from unit volume; the two terms a(Cco) and b(Co>) correspond respectively to the 
two source terms ('self noise' and 'shear noise') on the right hand side of Eq.(l 1). Cco is 
an 'effective' Doppler shifted source frequency 3,4 , and the outer factor C disappears in the 
integral over co (cf. also remarks under Eq. (52)). 

The format of Eq.(61) was effectively confirmed in the work of Pao and Lowson 24 . 

Their approach was, in fact, superior to Ribner's in that it yielded the convective 
amplification factor C" 3 in an internally consistent fashion. Ribner, on the other hand, 
used the 'compact eddy' assumption in a compatible M— >0 approximation to derive the 
Basic Directivity. Then he generalized this to finite M by assuming that the broadband 
factor C' 5 applied as a multiplier, to yield Eq.(61). 

The shear noise, given here by the b(Cco) term, has a dipole-like directivity. Thus it has 
been argued that the convective amplifier should be the C' 3 of a dipole, rather than the 
C' 5 of a quadrupole. Pao and Lowson 24 obtain a C' 3 multiplier, but they attribute the 
change from -3 to -5 to the choice of turbulence model. "The spectrum in the low wave 
number region follows a k 2 law, which artificially raises the power dependence from [-3 
to -5]. In actual measurements of jet turbulence, the spectrum in this region is usually 
flat." We observe, however, that the k 2 law refers to a 3D spectrum, whereas the hot 
wire measures only a ID spectrum, and the flat ID spectrum at small k is fully 
compatible with a k 2 3D spectrum. Thus there is no artificiality here, and the convective 
amplifier C -3 would appear to be applicable to the shear noise as well as to the self noise. 
Comparisons with experiment on this basis seem to support that conclusion. 

To our approximation, the magnitudes of C" 3 and IG/vl 2 apply equally to unit volume and 
to the jet as a whole. This is not true for the a(Cto) and b(Cto) of Eq. (61): the magnitudes 
differ for the two applications, but the ratio b/a is unaltered. An important prediction of 
the theory 4,23 that connects these two spectra is 

b(Cci)) = p a(2Cto) (62) 
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where (3 is a constant of order 2. The implications of this are elaborated below. 

The 'basic directivity’ of Eq.(61 ) is compounded of b(Cco) and a(Cco) in different 
proportions, depending on direction; credibility is afforded by comparisons with 
experiment. First, the predicted overall directional pattern over a range M=0.37 to M=0.9 
(outside the ’cone of silence’) is about right^.29 ; this is a minor point, since the 
contribution of the basic directivity is weak compared with that of the convective 
amplification. Second, and much more significant, the prediction in Eq.(62) of separate 
spectra of "shear noise' b(Cco) and 'self noise' a(2Cco) of almost identical shape, but with 
one octave relative shift, was strikingly confirmed^ 1 ). These are extracted from 
measurements at angles of 45 and 90 degrees from the jet axis, without adjustable 
constants. When normalized to unity peak, curves of shear noise b(C(0) and self noise 
a(2C(0) virtually collapse on one another. 

The third point has to do with correlations of two microphones located on a large sphere 
centered on the jet nozzle: with one microphone fixed, the other was displaced either 
along a meridian or a circle of latitude. It was with a series of such measurements carried 
out by Maestrello 41 that the theory was compared. In particular, correlations of two 
microphones are sensitive to details of the source instantaneous directivity, whereas the 
single microphone mean square response is not. The Lighthill-Ribner theory 23 leading to 
Eq.(61) was extended to deal with this case 42 ’ 4 -*. It was found that prediction of two- 
microphone cross-correlations along circles of latitude showed good qualitative 
agreement over a range of angular separations, and for different latitudes. This was true 
both in broadband 42 and the more demanding narrow band 43 ; see also references 51 and 
52. (Microphones located along a meridian exhibit a cusp-like correlation in broadband, 
decaying sharply with separation 41 . The agreement of the theory here was particularly 
striking 42 ; although only weakly relevant to Eq. (61), it makes a strong case for the 
credibility of key aspects of the basic theory 23 .) 

Basic Directivity x Convective Amplification x IG/vl 2 

The product of the weak basic directivity and the strong convective amplification, C* 3 , 
Eq.(61), yields the directivity of the jet noise of frequency 0) on the Lighthill-Ribner 
model. Multiplication by the normalized Green’s function, IGjyl 2 , modifies this to allow 
for the flow-acoustic interaction (refraction 3 * 3 and shielding 7 * 1 T 12 * 44 ). IG^yl 2 dominates 
in a cusp-like fashion near the axis to produce the refractive downstream ’cone of 
silence’. The broad fan-shaped amplification, C* 3 , with its maximum downstream as 
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well, amplifies not far from uniformly across the 'valley', scarcely modifying its shape. 
This accounts for the relatively close match between the measured IG^I 2 and the jet noise 
pattern near the axis (Fig.l). The resultant of the three factors is the well known heart- 
shaped directional pattern of jet noise. The opposing effects of C" 5 and IG/yl 2 , in their 
relative strengths, determine the location and magnitude of the maxima. The progression 
of patterns from Basic Directivity to Basic x C' 5 to Basic x C 5 x |G^| 2 is schematically 
shown in Figure 3; (the factors are additive on the decibel scale of the figure). 

Directivity vs Spectrum 

We have been discussing the directivity of jet noise in frequency bands; that is, the 
variation of <J>(xlo)) with d, the angle of the observer vector x with the jet axis. The 
directivity of the mean square sound pressure p^(x), being an integral over to, is a 
weighted average of these. The spectra at fixed angles i) are, of course, cross-plots. 

The directivity of the emission from a typical unit volume @ 3 d>(xlco)/3y 3 and its integral 
over co, 3 3 p 2 /dy 3 ) and the directivity of the entire jet, discussed above, are much the same 
(normalized to values at $ = 90 degrees). But the spectral shape emitted from unit 
volume differs greatly from that of the entire jet. The former is estimated as less than 2 
octaves wide, the latter is measured 54 as some 5 octaves, taken between the 6 dB down 
points. Thus calculations for a single unit volume may serve for the directivity of the 
entire jet, but an integration over the jet is required for estimating the spectrum. The 
directivity is only slightly sensitive to the asumed turbulence model (it affects the basic 
directivity); the spectrum requires a detailed estimate of turbulence properties throughout 
the jet. These refer to shapes : prediction of absolute levels requires a higher level of 
accuracy in the turbulence data and its modelling, and in any simplifying assumptions. 

APPROXIMATE QUANTITATIVE PREDICTION 
We turn now to the approximate quantitative prediction of jet noise properties. The 
central element is the four dimensional Fourier transform (4DFT) of the source term 
correlation function Rk m (e.g., Eqs. (3) and (53)). For illustration we limit attention to the 
self noise term, the shear term involving rather more complication. Neither theory nor 
measurements (so far) adequately describe the correlation function over space and time. 

If it is dealt with at all (in place of scaling arguments), a simplistic form (e.g., a 
Gaussian) is often assumed to ease the mathematics. It follows that the space transform 
(3DFT: »k) is of dubious accuracy, and the mathematical difficulty it entails seems 

thereby unwarranted. Thus the 3DFT is normally bypassed (an exception is ref. 24) by 
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means of the 'compact eddy' assumption: this reduces the transform to just a volume 
integral (over ^), as noted earlier. 


The 4DFT is then a volume integral times a 1DFT (x— ko). Moreover, the 1DFT defines a 
spectrum between one and two octaves wide. For simplicity this is commonly 
approximated as a single line: a 5-function of local frequency (Osc and amplitude A. 
Combining these yields the rather drastic approximation 

4DFT (correlation function) « A((uj c ^) 2 LiL 2 L 36 (o)-co S e)) (62) 

where the volume integral (of the normalized correlation function) is by definition the 
product LiL, 2 L 3 of the three orthogonal scales of turbulence; normalization is by the peak 
correlation (uk 2 ) 2 ), and Uk refers to the component of turbulence in the k-direction 


(aligned with the x-direction). All these parameters are functions of y. Insertion in 
Eq.(54) yields the power spectral density emitted from unit volume at y as 


j~ a 3 <E(xlci)) j 

8y3 


VM, SE : 


Po^Uk 5 ) 2 
‘ C 5 c 0 4 (4tix)2 


LiL2L35(G>-G)se) 


as applied to the self noise. This refers to a virtual medium at rest. The format of Eq.(53) 
generalizes this to a jet flow, encompassing refractive and shielding effects, via the factor 
\G N \ 2 (using Eq.(56)): 

[- ^(xl^ lsE” A\G^ 2 L iL 2 L 3 5(0)-tD se ) (64) 

dy i C 3 c 0 4 (4jtx) z 

Mathematically, the appropriate choices for the amplitude A and local frequency 0) S e are 
these: (i) A should be the integral of the 1 DFT spectrum (cf. above Eq.(62)) over to, and 
(ii) o)se should be the centroid or first moment of that spectrum. Then the 5-function 


approximation will yield minimum error in computing the overall power spectral density 
<t>(xlco) by integrating Eqs. (63) or (64) over y. In practice, heuristic scaling laws are 
often used for evaluating both A and (0 SC - These and the other parameters must be 
estimated as y ranges throughout the jet. A representative example 27 of such an 
estimation, leading to rather good spectral prediction, is given in Appendix D. 


CONCLUDING REMARKS 

Central Result 

The central result of the paper may be restated in simplified terms. Lighthill posed his 
aerodynamic sound sources as radiating into a ‘virtual medium at rest.’ Refraction of 
sound (creating the axial ‘cone of silence’) was suppressed by approximating the density 
in the dominant source term as constant (p = p 0 ). But by deferring the step p = p 0 , we 
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can pose the radiation as being emitted into the actual jet flow. This brings the refractive 
effect of the flow into play. Moreover, the residual sound source term is the same. 

Mathematically, the only change is replacement of the solution for a pure tone point 
source in a medium at rest by the solution for the source in a jet flow. The former can be 
written down by inspection as e lkr /47tr, the latter is a complicated solution, G(x,y|co), of a 
convected wave equation. But, at frequencies characteristic of jet noise, we find from 
both experiment 32 - 33 and calculation 8 - 9 that G reduces in the far field to e' kr /4rar (with a 

phase shift) times a directional factor. That directional factor, for a single round jet, is 
near unity for angle i3 greater than some value i3m- For smaller angles it decreases 
sharply to a minum on the jet axis, i) = 0. This describes the 'cone of silence' (Figs. 2,3). 
In summary, in the far field the new G differs in amplitude from the Lighthill e‘ kr /47tr 
significantly only within the ‘cone of silence’. Use ofG thus serves to extend the 
Lighthill-based solutions into this refractive zone. But outside it may he dispensed with, 
with little error. 

Relationship to Other Approaches 

A variety of Lighthill-based solutions — formalisms for jet noise prediction — have 
been used, e.g.. Refs. 23-27. They were all approximations. As discussed, they usually 
involved simplistic replacements for the four dimensional Fourier transform formalism. 
The turbulence correlation function, if it was modeled at all (rather than bypassed by 
heuristic assumptions), was normally taken as separable in space and time (implicit in 
Eqs.(62)-(64)). Chu 3() improved the model, removing the spurious separability: he used 
data from his own comprehensive program of very careful measurements by hot wire. He 
did, however, avoid the 4DFT by invoking the 'small eddy' assumption. Neverthless, 
despite these deemed improvements, his predictive accuracy fell far short of the best of 
refs. 23-27. Hindsight suggests the capabilities of his data may not have been optimally 
exploited. A revisit in the light of the present formalism could be profitable. 

The predictive problem is compounded by the difficulty of a four-dimensional Fourier 
transform. As noted, the 'compact eddy' approximation reduces this to a simple volume 
integral. But that step compromises the accuracy of the prediction of convective 
amplification. The approach has, indeed, led to fairly accurate predictions, both for round 
and more complex jet configuration. However, even with this approximation, they could 
be improved as indicated under Eq.(64). The further inclusion of the factor IGyvl 2 will 
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extend the solution into the 'cone of silence’ (small 0 region), and even improve the 
accuracy outside this region (by the amount IG/vl 2 differs from unity). 

Lilley’s wave equation, in the hands of Balsa, 15 leads to a similar formula IG #| 2 x four- 
dimensional Fourier transform of source correlation function. Here again, in practice the 
Fourier transform is bypassed by an approximation. There are other important 
differences . A major component, the shear term included in Eq.(l 1), is missing from the 
source term: the consequences are discussed in the text. And the squared normalized 
Green’s function, IG#| 2 , refers to a moving point source. The effects of source convection 
and wave convection — respectively governing amplification and refraction (via velocity 
gradients) — are thereby combined. In contradistinction, these effects are decoupled 
herein by the use of the value of IGyyl 2 for a stationary source That is, IG #| 2 plays 
essentially the same role as C" 5 IG/vl 2 : they should be largely equivalent. 

It is this decoupling that allows the simplicity of the Lighthill-based formalisms to be 
applied outside the ‘cone of silence', since IG/vl 2 is near unity there. Another advantage is 
that this refractive IGyvl 2 is determined for a realistic spreading jet — either by 
calculation or experimentally — whereas IG #| 2 has been evaluated only for idealized, 
infinite, nonspreading jets. 

Issue of Shielding 

Mani 11 * 12 anc i Balsa 15 have pointed to a 'shielding' role of the mean flow in reducing the 
convective amplification (a function of direction) at high frequencies. But 
comparisons 28 * 29 with the directivity measured by several investigators do not bear this 
out. The simple theoretical convective factor C ' 5 adequately predicts the curves up to 
M=0.9, source Strouhal no. 1.0 (observed Strouhal no. 1.33). This is a reduction in 
convective amplification compared with the classical factor (1- M c cos0)' 5 , but it is not an 
effect of flow shielding. Instead, it results from retarded time variation across an 'eddy', 
not allowed for in the classical factor. 

It would seem that the frequencies of jet noise are simply too low: for significant 
shielding the flow dimensions must be much larger than a typical wave length of the 
sound 55 * 5 ^. This is a requirement for geometric acoustics (ray acoustics) to apply. 
Schubert's calculations 8 * 9 show that jet noise is very far from that regime. 
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These remarks refer to a single round jet. The case of multitube jets (or equivalent 
corrugated nozzle jets) is another matter. Substantial shielding of the high frequency 
noise of the inner jets by a ring of the outermost jets is a demonstrated fact. Balsa 44 , via 
a moving source Green's function, shows apparent agreement with measurements that he 
cites. No stationary source Green's function, as proposed herein, has as yet been 
evaluated for this scenario. 

Range of Applicability 

The results herein are for the far field only. Moreover, they are presumed to beapplicable 
primarily for subsonic jets. For supersonic jets additional noise sources come into play. 
Tam 60 , in his review article, develops the case for instability waves, identified as 'large 
scale coherent structures', being a major source of noise. He evaluates the noise directly 
via a 'stochastic wave model' with very impressive agreement with experiment. But we 
note that these instability waves, to the extent that they coexist with the random 
turbulence, will contribute to the correlation function Rk used herein. But evaluation is 
another matter: the similarity laws for jet turbulence, which have been used with success 
for subsonic jets, would have to be reevaluated for the supersonic regime. 

Perhaps more importantly, there will be a pattern of shock waves if the jet does not issue 
at the design speed from a properly contoured convergent-divergent nozzle. It was shown 
many years ago that shock-turbulence interaction would generate intense noise^ 7 '^. j n 
recent years Tam (summarized in ref. 60) has attempted quantitative prediction of this 
shock-associated noise; he analyzed the interaction between instability waves and a 'wave 
guide' model for the shock structure. His near field patterns show a close match to 
measurements. 

The results are further restricted to jets issuing into ambient fluid at rest: that is, static test 
conditions. The effects of forward flight on the jet noise are not considered. Michalke 
and Michel29.49.50 have extended the Lighthill theory to provide a successful prediction 
of these effects. This takes the form of a scaling law that maps the intensity of a static jet 
at certain jet Mach number and direction into that for a moving jet at an altered Mach 
number, direction, and distance. (See also an approach via CFD methods 61 .) Refraction, 
governing extension into the 'cone of silence', is not allowed for: this would involve a 
further development of the present stationary-source Green's function approach. 
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APPENDIX A: GENERALIZED CONVECTIVE WAVE MODIFICATIONS OF 

LIGHTHILL EQUATION 


EXACT WAVE EQUATION 
Expansion of Lighthill Source Term, Ql 

The restriction to a transversely sheared mean flow of uniform density, Eq. (8), is relaxed 
here: Lighthill's source term expression, which we shall call Ql (the right-hand side of 
Eq. (5)), is expanded under the specifications 

v i = Uj + uj; Ui = Ui(x); (vi)av = Ui (Al) 

p = p(xj) + p’ ; (p) av = P(xj) (A2) 

By Csanady's 7 Eq. (3), Ql expands as follows: 


a2pv 1 v j+ _i_a2 E ^av i+ ay i 9y i \ + _L^E 

5 xi 0 xj + c 0 2 dt 2 dt 2 p [ 5 xj 5 xj + 5 xi 5 xjJ + Co 2 3 t 2 


02p a2p 92p 

a t 2 V ‘ c)xic>t ' V|V j dxjdxj 


(A3) 


With this expansion it is easily shown that Eq. (5) of the main text is equivalent to 
Schubert's ^- 9 exact wave equation for an inviscid nonheatconducting fluid. Inserting 
Eqs. (Al) into (A3) yields 7 

„ /aujaui auiau;\ a 2 p OTT a 2 p TTTT a 2 p 
p (a Xj a Xj + axj axj J J’ at 2 " ‘ 5xj5t" ' j axjaxj 


i a 2 p a 2 pujUj aujapuj a r au, ~\ 

+ c 0 2 dt 2 + dxjdxj + 2 a Xj a Xi + 2 a XJ [ pu i axj J 

When Uj = (U(x 2 ), 0, 0) this reduces to the right-hand side of Eq. (8) of the main text. 
With the further expansion 

a 2 p Ui uj , aujapui a 2 uiui » auiam « aujap .aujujap a 2 p 

axjdxj + a Xj axj - ^axjaxj + ^axjaxj + axj 5x[ + axj axj + u ' u j ax^axj 

(A5) 

and the definition 

£? a 2 a 2 a 2 

^ s § + 2u ^ + u i u i 5 ^ < A6 > 


and some rearrangement, Ql becomes 
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Ql = 2p 


3Ui 3uj 

3xj 3xj 


1 3 2 p t^ujuj _ 3Uj 3p 

c 0 2 5t 2 Dt 2 3xj5xj + u -» 9xj 8xi 


9ujUj 3p 
+ z 3xj a Xi 


0 2 p 

+ UiU Ja^j 


2 — 

2 a Xi 





/a^aui + au]auj \ 
l axj a Xj a Xi ax| I 


(A7) 


This expansion of the acoustic source terms, the right-hand side of the wave equation (5), 
is exact. 


APPROXIMATE WAVE EQUATIONS 
Incompressible Turbulence 

We now introduce approximations in two stages. First, in all but the first three terms of 
Ql, we neglect the compressibility of the turbulence in the application to subsonic flows. 
We argue that density perturbations p' in the nonexcepted terms account for scattering of 
sound by turbulence, and they may be neglected in dealing with generation of sound. 
Thus, in these terms, p may he replaced by its local temporal mean p(x). Consistently, it 
is implied that uj, in all but the excepted terms, contains no compressible component. 
Despite this assumed incompressibility of the source terms, sound (pressure waves) will 
indeed be generated as Lighthill 1 ’ 2 showed. Thus compressibility has been retained 
where acoustically necessary: in the final left-hand-side wave operator. 

The first excepted term is 2p(dUi/9xj)(duj/c)xj). By the argument in the main text above 
Eq. (1 1), the very small compressible, or acoustic part of this term, being linear in Uj, 
participates to the first order in wave propagation; thus we move it to the left-hand-side 
wave operator. The terms c Q 2 d 2 p/dt 2 and -6 2 p/Dt 2 have also to do with wave 

propagation; we move them likewise to the left-hand side, approximating ft^p/Dt 2 by 
C' 2 6 2 p/Dt 2 (Appendix B). 

With these term shifts from the right-hand side (Ql) and approximation of p by p therein, 
Eq. (5) may be rewritten. First we need an approximate equation for mean flow 
continuity. Taking the correlation p'u; as negligible yields this as 

dpUj/axj = 0 (A8) 

By virtue of Eq. (A8), two terms of the approximate Ql may be collapsed into one: 
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The modified Eq. (5) then reads 
1 


n 


l^E + 


3 2 p . nit 9 2 p , dUj 3uj ac n'y -d 2 UjUj . „-dUi 

2Ui 5Jj3t + U ' U j3x^ 1 - -mt - 72 P - Pdxjifej + 2 P^ 


j 3uj 


+ 2 


dujUj 3p 


j 3xj dxj 9xi 


3 2 p 

+ UiU J5^rj 



(A10) 


Transversely Sheared Flow 

Let us specialize now to a transversely sheared flow (Eq. (7) of the main text), with 
transverse density gradient as well: 


U i = (U(x 2 ), 0, 0); p = p(x 2 ) 
The modification (A10) of Eq. (5) simplifies to 


(All) 


c 2 


.au2Ui ap 


+ U2 2 


3 2 p 

5x2 2 


But it can be quickly verified by direct expansion that 

- 3 2 VjVi - 3 2 UjUj 

P 3xi3xj ~ P 3xj3xj + 
so that an alternative form is 


^3U3u 2 

2p 5x25xT 


(A 12) 
(A 13) 


J_r^R 

— 2 3t 2 

C l 


+ 2U 


a 2 p 

dxi3t 




-au au^ 

' zp ax 2 axi 



3u2Uj dp 
3xj 3x2 


+ U2 2 


a 2 p 

3x2 2 


(A 14) 

where it is recalled that vj is the instantaneous resultant flow Uj + uj, as defined in Eq. 
(Al). An equation of the form (A 14) results also for a cylindrical jet Uj = (U(r), 0,0), p = 

p (r), if U2, X 2 are replaced by u r , r, respectively (see Eq.(A15). Equations (A12) and 
(A14) are the key results of this Appendix; they generalize Eqs.(ll) and (13) of the main 
text, respectively, to the case of flows of nonuniform mean density. The applications and 
implications are discussed below. 


Jet Flow 

Following Schubert 8 - 9 we note that Uj in a jet is essentially unidirectional: 



Ui = [U(r,xi), 0, oi; r = Vx 2 2 + x 3 2 round jet (A15) 

= x 2 two dimensional jet 

with a strong dependence on r, and a weak dependence on xj. Further, 

p = p(r;xj) (A 16) 

with a similar dependence. Thus, for both U and p, the gradients along \\ are very much 
less than those along r. For the foregoing transversely sheared flow the xi-gradients are 
identically zero; requiring this led to Eqs. (A 1 2) to (A14). It follows that Eqs. (A 12) to 
(A 14) may be applied to jet flow as a close approximation, with x 2 replaced by r. 


Density Scenarios: p = constant vs. p = (5 (r) 

Let us specialize further to a uniform mean density, 

p = p 0 = constant (A 1 7) 

This, together with Eq. (All), recovers the scenario of the main text. It is seen that the 
density gradient terms drop out, and Eqs. (A 1 2) and (A14) reduce to Eqs. (11) and (13), 
respectively: the modified Lighthill equation in the form ofEq.(13 ) is confirmed as a 
special case of the more general form of this Appendix, Eq.(A14). 


From the foregoing, it is clear that density gradients, via the additional source terms, 
cause more noise to be generated. This has been explored in the context of hot jets by 
Morfey 2 ^, Mani 12 , and by Michalke and MicheH^O. The source terms in (A12) and 


(A 14) appear similar to those deduced by Mani. The term 2 dg 2U j is essentially of 

dipole form, 3Qj/3xj, (treating 3p/3x 2 as a spatial constant): it would yield a factor ltd 2 in 
place of IkI 4 in an equation like (4). As a consequence of this, or by arguments given in 
the cited references, the corresponding radiated sound power would vary as U 6 ; they 
showed it could exceed the ordinary quadrupole-source jet noise, with its U 8 law, for 

sufficiently hot jets. The term u 2 2 3 2 p/3x 2 2 is of monopole form, leading to a U 4 law 
(ltd 4 factor replaced by unity). This would radiate very weakly, the curvature 9 2 p/3x 2 2 


being minimal in the zone of strongest turbulence, where the mean flow shear and 3p/3x 2 
both maximize. 


Over time "...dozens of equivalent (and nonequivalent) source term expansions have been 
published by flow noise researchers. This multiplicity of competing source terms has 
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been a major contributor to confusion..." (Reference 29). In this author’s view, the 
effective limit in the law of diminishing returns has been reached in the expansions of this 
Appendix: of the Lighthill wave operator (left hand side of Eq.(AlO)) and of the Lighthill 
source term for an inviscid nonheatconducting fluid (exact, Eq.(A7); approximate, 
Eq.(A14)). 


APPENDIX B: APPROXIMATION OF ft 2 p/Dt 2 as C' 2 5 2 p/Dt 2 


The exact relation for a moving fluid element of constant entropy 




Dp _ J_ Dp 
Dt q2 Dt 

(Bl) 

implies 


D 2 p 1 D 2 p Y-l /Dp>2 
Dt 2 c 2 Dt 2 p c 4 J 

(B2) 

where D 2 /Dt 2 is the convective second derivative 

D 2 d 2 , ,, d 2 d 2 , Dvj d 

Dt 2 dt 2 V ' dxjdt + V ' V j dxjdxj + Dt dxj 

(B3) 

in which 


Vi = Uj + Uj 

(B4) 


is the mean velocity Uj plus a perturbation uj. This may be compared with a mean flow 
convective derivative 

6 2 a 2 a 2 a 2 

5^ = P + 2Ui S5 + UiU i3^xj < B5 > 


It is argued, following Schubert 8 - 9 , that replacement of (B3) by (B5) (and c 2 by c 2 ) 
merely supresses scattering of sound by the jet turbulence; (Dvj/Dt>a/axj is neglected 
also, as yielding second order terms. 

We illustrate this in terms of a plane sound wave in a reference frame following the mean 
flow: 

p = p + Ip'l e‘( k i x i - wt); Vkj 2 + 1<2 2 + k 3 2 = k = oVc (B6) 

p = p + Ip'l e'( k j x i - “0 

(Herein the contribution of turbulent perturbations to c 2 , c 2 , andl/c * 2 are considered 
negligible, so the designated local time averages are taken to be equivalent. With this 
assumption , the averages involving soundspeed c in Eq.(B6) and later equations are 
compatible.) Using the momentum equation. 
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(B7) 


1 Dvifo _ 1 ,dp 

Dt 3xj " _,/dxj ' 
c z pc z 


And the amplitude is (with p 


P) 

1 Dvj3p 
Dt 3x; 

r l 


Jil|p.|2,i^ p .|2 


pc" 


TP 


(B8) 


For comparison we evaluate c* 2 D 2 p/Dt 2 , which is c' 2 3 2 p/3t 2 in the moving frame, to a 
close approximation. By (B6) the amplitude is 

|-L^U^|p’| = k 2 lp , | (B9) 

c 2Dt c 2 

The term (B8), arising from the operator (Dvj/Dt)3/3xj, is seen to have an amplitude a 
factor lp'1/yp smaller than the term (B9). This justifies neglect of that operator. 


By a similar example we can show that the last term of Eq. (B2) is of higher order and 
may be dropped. Thus, on the basis of these order-of-magnitude estimates, we conclude 


that 

6 2 p __ i 6 2 p 

Dt 2 ~ - 2 Dt 2 
c z 

is a sufficient approximation for the modified Lighthill wave equation. 


(BIO) 


APPENDIX C: APPROXIMATION OF K AS k 
For frequencies characteristic of jet noise, we examine the change in phase \}t in the far 
field when the source is displaced from y to y'. Equation (25) gives a linear 
approximation as 

|v|/' - Yl = [y(x,y'lco) - y(x,ylG))] = -K«(y’-y) (Cl) 

where 

K = — (Vy V|/)x»y (C2) 

and wave convection by the jet flow is allowed for. In the absence of flow (or for a 
virtual medium at rest) this reduces to 

k = - 1 V y (kr) 1 x »y = | kr/r] x » y » kx/x (C3 ) 

with corresponding phase change 

lV-VlvM = -k*(y’-y) (C4) 

This is an approximation to the exact relation 

W - y]v M = k(lx-y'l - Ix-yl) (C5) 
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The (negative) phase advance due to the flow is the left hand side of Eq. (Cl) minus that 
of Eq. (C4); designating this as S, it is obtained as 

S s W - V] - W - VlvM = (k-K)*(y'-y) (C6) 

Schubert 8 , in his numerical computations for a point source in a jet flow, has evaluated 
this phase difference in a number of cases (he used the more precise relation (C5) in 
defining the left hand side as 

S a W - yl - kflx-y'l - Ix-yl) (C7) 

effectively: for the large x/y ratio of our examples the difference is negligible). In his 
Figures 30 to 32 he plots S/7t vs (y’i- yi)/D; this corresponds to y' and y both directed 
along the jet axis; x is directed at angles fr= 0, 6.8, and 16.2 degrees. 


We take the simplest case, fr=0. Now k is directed along x by Eq.(C3), and when x lies 
along the axis, so does k. By virtue of the flow symmetry about fr=0, k must also be 
directed along the axis. For this scenario Eq.(C6) becomes 

S = (k-K)(y’i-yi) (C8) 

so that 


K/k = 1- S/k(y‘ i -y i ) = 1- e 

The ratio K/k is evaluated in several examples with data from Figs.30 and 32. The 
specifications are: Mj = 0.3, , lx-yl/D=100. The results are: 

W 0 ^kD= 1.055 and (y'i-yi)/D = 1: K/k = 1- jt(.005)/[ 1.055(1)1 = 1 -.015 

(y’l-yi)/D = 2: K/k = 1- 7t(.015)/[l .055(2)1 = 1 -.022 

Wo=kD= 2.46 and (yVyi)/D = 0.5: K/k = 1- 7t(.032)/(2.46(0.5)] = 1 -.08 

(y'l-yi)/D= 1: K/k =1- 7t(. 101 )/[2.46(l )] = 1 -.13 

It is seen that K/k differs little from unity: the difference le| increases with kD 


(C9) 


(proportional to frequency) and with (y'i-yi)/D. The values kD = 1.055 and 2.46 
correspond to the upper part of the jet noise spectrum, which has a broad peak around 
kD= 0.3. Therefore, for the bulk of the spectrum the difference le| will be less for a given 

(y'l-yi)/D. Furthermore, an appropriate source displacement (y’l-yi) must be less than 
the local macroscale L of the jet turbulence, a function of yj. Estimates are L/D = .35, 

.15 for kD = 1.055, 2.46, respectively. Thus the values of (y’ryi)/D in the examples are 
unacceptably high, especially at the highest kD. It follows that appropriate values of le| 
are substantially smaller than those of the examples: K/k differs from unity by not more 
than several percent at jet noise frequencies for Mj = 0.3. The difference will, of course, 
increase with jet Mach number. 
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On going to fr > 0, the curves of S/7t vs (y'i*yi)/D in Figs. 30 and 32 are not greatly 
different: the phase difference S between the flow and no flow cases remains very small 
as f) varies. This implies that the phase gradient (V y \|0 x »y for the flow and no flow 
cases has nearly the same direction as well as magnitude. That is, the inference from the 
d = 0 examples above, with respect to magnitudes, that 

k* k (CIO) 

may be generalized to the directions as well; that is, 

K = k (Cl 1) 

to a close approximation for at least low speed jets at typical frequencies. 


APPENDIX D: MOON-ZELAZNY MODEL OF THE TURBULENCE 

PARAMETERS 

Moon and Zelazny 27 derived an equation roughly equivalent to the approximate Eq.(63) 
for the self noise. (They had another equation, more complex, for the shear noise.) These 
were based on a space-time (rather than wave number- frequency) domain solution of the 
Lighthill equation and involved various simplifying assumptions. In particular, the 
operation 3/d t was interpreted as effectively multiplying by a characteristic frequency (a 
function of axial position yi in the jet): 0) sc for the self noise and ( l/2)co se for the shear 
noise. There was also an implicit heuristic assumption for the constant A of Eqs. (62)- 
(64). For implementation for the scenario of a round jet, they developed a model for 
determination of these and other needed turbulence parameters. A major feature (taken 
from a Ph.D. dissertation underlying ref. 27) was an eddy viscosity formalism for 
evaluation of the rms axial component of turbulence, u’, as a function of r and y\ 
throughout the jet. This led to evaluation of (u^ 2 ) 2 appearing in Eq.(63). Other heuristic 
features of their model are: 

L s = shear layer half-width L] = 0.358L S L2 = L3 = 0.179 L s 

C0s e = 0.3U(yi)/L2 o) sh = 0.3U(yi)/Li 

The agreement with both turbulence and noise measurements was good. These were 
absolute levels: in particular, for the noise there was no shifting of the decibel scale for 
best match . This has to be impressive. On the other hand, the assumptions, e.g., of 
scales and frequencies, involves empiricism. 
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Fig.1 . Shear term in source role. Match of experimentally derived shear noise and self 
noise spectra, normalized to same amplitude. The self noise has been downshifted one 
octave, in accordance with the theory. (After ref. 29, based on ref. 40.) 
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Fig. 2. Match of noise intensity patterns to indicate refraction-dominated zone. Left-hand 
side: room temperature jet (M=0.5) with reduced intensity ’cone of silence’ (rays turn 
outward); right-hand side; very cold jet (-180°C) with enhanced intensity lobe (rays turn 
inward to a quasi-focus). (After ref. 29, based on refs. 32 and 33). 


42 




B: Basic Pattern 
C: Convective 
Amplification 
R: Refraction 
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Fig.3. Convective amplification (C- 5 ) and refraction (normalized Green's function IGnI 2 ) 
modify the basic pattern of intensity vs direction (in decfoels). (After ref. 29, based on 
refs. 32 and 33). 
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